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We report the observation of the decay mode B — > (jxf)K based on an analysis of 78 fb _1 of data 
collected with the Belle detector at KEKB. This is the first example of a b — > slsss transition. The 
branching fraction for this decay is measured to be B(B ± — » <j)(j)K ± ) — (2.6^0 9 ± 0.3) x 10~ 6 for a 
4>4> invariant mass below 2.85 GeV/c 2 . Results for other related charmonium decay modes are also 
reported. 

PACS numbers: 13.25.Hw, 14.40.Nd 



We report the observation of the decay mode B — > 
4>4>K, the first example of a b — > sssss transition. In 
the Standard Model (SM), this decay channel requires 
the creation of an additional final ss quark pair than in 
b — > sss processes, which have been previously observed 
in modes such as B — > ipK . In addition to improving our 
understanding of charmless B decays, the (jxfi K state may 
be sensitive to glueball production in B decays, where 
the glueball decays to c/xf> . In addition, with sufficient 
statistics, the decay B — > (fxf)K could be used to search 
for a possible non-SM CP-violating phase in the b — > s 
transition |2| . Direct CP violation could be enhanced to 
as high as the 40% level if there is sizable interference 
between transitions due to non-SM physics and decays 
via the rj c resonance. 

We use a 78 fb _1 data sample collected with the Belle 
detector at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider [3j operating at the T(45) resonance 
(Vs = 10_58 GeV). The sample contains 85.0 x 10 6 pro- 
duced BB pairs. The Belle detector is a large-solid-angle 
magnetic spectrometer consisting of a three-layer silicon 
vertex detector, a 50-layer central drift chamber (CDC), 
a system of aerogel threshold Cerenkov counters (ACC), 
time-of- flight scintillation counters (TOF), and an ar- 
ray of CsI(Tl) crystals located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to 
identify K Q L and muons. The detector is described in 
detail elsewhere 0. 

We select well measured charged tracks that have im- 



pact parameters with respect to the nominal interaction 
point (IP) that are less than 0.2 cm in the radial direc- 
tion and less than 2 cm along the beam direction {£). 
Each track is identified as a kaon or a pion according to 
a K/ir likelihood ratio, Ck/{£tt + Ak')j where CkM are 
likelihoods derived from responses of the TOF and ACC 
systems and dE/dx measurements in the CDC. We se- 
lect kaon candidates by requiring £k/(Ctt + C-k) > 0.6. 
This requirement is 88% efficient for kaons with a 8.5% 
misidentification rate for pions. Kaon candidates that 
are electron-like according to the information recorded 
in the CsI(Tl) calorimeter are rejected. 

Candidate <f> mesons are reconstructed via the <f> — > 
K + K~ decay mode; we require the K + K~ invariant 
mass to be within ±20 MeV/c 2 (±4.5 times the full 
width) of the mass 0. For the B°(B°) -> (j)<j>K% de- 
cay mode, we use Kg — ■» ir + ir~ candidates in the mass 
window 482MeV/c 2 < M{tt+tt-) < 514MeV/c 2 (±4cr), 
where the distance of closest approach between the two 
daughter tracks is less than 2.4 cm, the magnitude of the 
impact parameter of each track in the radial direction ex- 
ceeds 0.02 cm, and the flight length is greater than 0.22 
cm. The difference in the angle between the pion-pair 
vertex direction from the IP and its reconstructed flight 
direction in the x — y plane is required to be less than 
0.03 radians. 

To isolate the signal, we form the beam-constrained 

mass, M bc = \Z^boam ~ Brecon | 2 , and the energy dif- 
ference AE = S rec on ~ ^bcam- Here £ b cam is the 
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FIG. 1: (a) <fxj> invariant mass spectrum. The open his- 
togram corresponds to events from the B — > (p(f>K ± signal 
region and the shaded histogram corresponds to events from 
the Mb c -AE sidebands, (b) M K + K - of one cf> meson can- 
didate versus M K + K - of the other for the events satisfying 
< 2.85 GeV/c 2 . The events are concentrated near the 
point (M^jMc/y). The dashed box shows the region selected 
for the B — » <j><j>K analysis. 



beam energy, and E Te 
energy and momentum of the signal candidate in the 
T(45) center-of-mass frame. The signal region for AE 
is ±30 MeV which corresponds to ±3.1<r, where a is 
the resolution determined from a Gaussian fit to the 
Monte Carlo (MC) simulation. The signal region for 
Af bc is 5.27GeV/c 2 < M bc < 5.29GeV/c 2 . The beam- 
constrained mass resolution is 2.8 MeV/c 2 , which is 
mostly due to the beam energy spread of KEKB. 

The major background for the B — > (fxpK process is 
from continuum e + e~ — > qq production, where q is a light 
quark [u, d, s, or c). Several event topology variables are 
used to discriminate the continuum background, which 
tends to be collimated along the original quark direction, 
from the more isotropic BB events. Five modified Fox- 
Wolfram moments, the S± variable and the cosine of 
the thrust angle are combined into a Fisher discriminant 
@- We form signal and background probability density 
functions (PDFs) for this Fisher discriminant and for the 
cosine of the B decay angle with respect to the z axis 
(cos#b) for the signal MC and sideband (5.20GeV/c 2 < 
A/ bc < 5.26GeV/c 2 and 0.1 < \AE\ < 0.2 GeV) data, 
respectively. The PDFs are multiplied together to form 
signal and background likelihoods, £ s and £ BG . The 
likelihood ratio £1Z = £s/(£s + £bg) is then required to 
be greater than 0.1. This requirement retains 97% of the 
signal while removing 55% of the continuum background. 

Figure ^a) shows the <jxj) invariant mass spectrum for 
events in the B^ — > (fxfrK signal region, where a clear 
r] c peak and some excess in the lower mass region are 
evident. 

To extract signal yields, we apply an unbinned, ex- 
tended maximum likelihood (ML) fit to the events with 
\AE\ < 0.2 GeV and M hc > 5.2 GeV/c 2 . The ex- 
tended likelihood for a sample of N events is £ = 



-(Ns+N B )Y\ N 



where vf" B ^ describes 



the probability for candidate event i to belong to the sig- 
nal (background), based on its measured M bc and AE 
values. The exponential factor in the likelihood accounts 
for Poisson fluctuations in the total number of observed 
events N. The signal yield N$ and the number of back- 
ground events Nb are obtained by maximizing £. The 
statistical errors correspond to unit changes in the quan- 
tity x 2 = — 21n£ around its minimum value. The sig- 
nificance of the signal is defined as the square root of 
the change in x 2 when constraining the number of signal 
events to zero in the likelihood fit; it reflects the prob- 
ability for the background to fluctuate to the observed 
event yield. 

The probability V for a given event i is calculated as 
the product of independent PDFs for M bc and AE. The 
signal PDFs are represented by a Gaussian for M bc and 
a double Gaussian for AE. The background PDF for 
AE is a linear function; for the M bc background we use 
a phase-space-like function with an empirical shape ||. 
The parameters of the PDFs are determined from high- 
statistics MC samples for the signal and sideband data 
for the background. 

For M(<fi(f)) < 2.85 GeV/c 2 , the region below the charm 
threshold, the ML fit gives an event yield of 7.3^2' 2 
with a significance of 5.1 standard deviations (a). Pro- 
jections of the AE distribution (with 5.27GeV/c 2 < 
M hc < 5.29 GeV/c 2 ) and of the M hc distribution (with 
| AS | < 30 MeV) are shown in Figs. EJa,b). As a consis- 
tency check, a ML fit to the projected AE distribution 
(Fig- EJb) only gives a signal yield of 7. 5^2 ^ with a 4.8cr 
statistical significance. Figure^b) shows a scatter plot of 
the two K + K~ invariant masses for events in the B me- 
son signal region with M(K + K~K + K~) < 2.85 GeV/c 2 
with the <f> mass requirements relaxed. Here there is a 
clear concentration in the overlap region of the two <j) 
bands. There is no event excess in the 4> mass sidebands, 
which leads us to conclude that the observed signal is 
entirely due to B ± — > (f><pK ± . Using a signal efficiency 
of 3.3%, obtained from a large-statistics MC that uses 
three-body phase space to model the B^ — > <j><j>K^ de- 
cays, we determine the branching fraction for charmless 
B ± -> 0(/.S ± with M u < 2.85 GeV/c 2 to be 

B{B ± -> #i?f ± ) = (2.6±J;J ± 0.3) x 10~ 6 , 

where the first error is statistical and the second is sys- 
tematic. 

Contributions to the systematic error include the un- 
certainties due to the tracking efficiency (5.4%), parti- 
cle identification efficiency (5%), and the modeling of 
the likelihood ratio cut (2%). The error due to the 
modeling of the likelihood ratio cut is determined us- 
ing B~ — * D a (-^ K~ir + Tr~ir + )ir~ events in the same 
data sample; these events have the same number of final- 
state particles and an event topology that is similar to 
the B ± — > (fxjiK^ 1 signal. The uncertainty due to the MC 
modeling (4%) accounts for the dependence of 
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FIG. 2: Projections of Mbc and A_B overlaid with the fitted 
curves for (a, b) B ± -> ^A"* with Af 00 < 2.85 GeV/c 2 , (c, 
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(e, f) _B -> rjcif* and r? c 



2{K + K~), and (g, h) _B -> J/^lf* and J/V> -> 2(K + K~). 



the detection efficiency. The systematic error in the sig- 
nal yield (6%) is determined by varying the means and a 
of the signal and the shape parameters of the background. 
We determine an upper limit of 5% on the possible con- 
tamination by non-resonant i? ± — > <j>{K + K~)^K ± or 
B ± — > 2(K + K~)nrK ± decays by redoing the fits with 
the <fr mass requirement relaxed. The sources of system- 
atic error are combined in quadrature to obtain the final 
systematic error of 12%. 

For the B°(B°) — > (f>4>Kg mode, there are only four 
signal candidates. We combine the S ± — > cjxpK^ and 
B°(B°) — ► 4>4>Kg modes and perform a ML fit and ob- 
tain a signal event yield of 8.71?,' g with 5.3<r statistical 
significance. Assuming isospin symmetry, we obtain 

B(B -> 4>4>K) = (2.3t{|J ± 0.3) x 1(T 6 , 

for M H < 2.85GeV/c 2 . 

No enhancement is observed in the region corre- 
sponding to the f j (2220) glueball candidate |5j, also ref- 
ered to as £. Assuming the mass and width of /,/(2220) to 
be 2230 MeV/c 2 and 20 MeV/c 2 , we define a signal region 



of 2.19GeV/c 2 < M H < 2.27GeV/c 2 , 5.27GeV/c 2 < 
M bc < 5.29GeV/c 2 and \AE\ < 30 MeV. One event is 
observed in this region with an expected background, es- 
timated from the sideband, of 0.5. Using an extended 
Cousins-Highland method that uses the the Feldman- 
Cousins ordering scheme and takes systematic uncertain- 
ties into account 0, we obtain a 90% confidence level 
(CL) upper limit of 3.7 signal events, which corresponds 
to 



B(B ± -> f J (2220)K ± ) xB(f >(2220) 



< 1.2x10" 



We select B 



± 



i] c K ± , r\ c 



quiring 2.94GeV/c 2 < < 



> 4>(j) candidates by re- 
3.02GeV/c 2 . A clear 



signal is evident in Figures |2fc,d), and the fitted yield 
of Ns = 7.0^2 3 events has a significance of 8.8cr. The 
corresponding branching fraction is 



B(B ± -> r ?c X ± ) x B{ Vc - 

= (2.2to° ± 0.5) x 10~ 6 . 

In addition to the previously listed sources of systematic 
error, here the error also includes the possible contam- 
ination from charmless B^ — > cfxj)K ± decays, which is 
estimated to be less than 1.2 events. Using the measured 
branching fraction B(B ± -> r/cK^ = (1.25 ± 0.42) x 
10 -3 fiof . we determine the r\ c — > <f>4> branching fraction 
to be 



l.8±o.t ± °- 7 ) x 10 ~ 



which is smaller than the current world average value of 
(7.1 ±2.8) x 10" 3 0. 

Since the J/ip and rj c charmonium resonances also de- 
cay to 2(K + K~), the decay chains B — * charmonium+J-T 
with charmonium — > 2(K + K~) can provide consistency 
checks of the B — > (jxfiK analysis. To select B — » 
2(K + K~)K candidates, we apply tighter particle iden- 
tification and continuum suppression requirements than 
in the case of B — > (fxpK in order to reduce the larger 
combinatoric background. Figure 01a) shows the invari- 
ant mass distribution of any two pairs of K + K ~, M^k, 
between 2.8 GeV/c 2 and 3.2 GeV/c 2 for the events in the 
B signal region. Significant contributions from both r\ c 
and J ftp intermediate states are seen. 

To identify the signals from rj c and J/tp interme- 
diate states, we require that the invariant mass of 
2(K+K~) satisfy 2.94GeV/c 2 < M iK < 3.02GeV/c 2 
and 3.06GeV/c 2 < M iK < 3.14GeV/c 2 , respectively. 
We use signal yields from ML fits to determine branching 
fractions. Figures Ele-h) show the M^ c and AE projec- 
tion plots with the fitted curves superimposed. Table [I] 
summarizes the signal yields, efficiencies, statistical sig- 
nificances, and the branching-fraction products. By re- 
quiring the invariant mass of one of the K + K~ pairs to 
correspond to a 4> meson, we also measure the decays of 
B ± -> T] c {J/ij;)K ± and r) c (J/ip) -> 4>K+R-. The results 
are included in Table [I] 
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TABLE I: Signal yields, efficiencies including secondary branching fractions, statistical significances and branching fractions of 
B — > (f><f)K and related decays. The branching fractions for modes with K + K~ pairs include contributions from — * K + K~ . 

Mode Yield Efficiency (%) Significance (a) B (xlO -6 ) 

B ± -> (M^ < 2.85 GeV/c 2 ) 7.3 

B -> 00ff (M 00 < 2.85 GeV/c 2 ) 8.7 ±!j;g 
B* -> ,/j(2220)A' ± , /j(2220) -> 00 < 3.7 

B^ — > r/cK^ 1 , r) c — > cf>(j> 7.0 ± 3 2 ° 3 

B ± -^r] c K ± ,r] c ^<j}K + K- 14.1 t 4 .^ 

B ± ^ Vc K ± ,r lc -^2(K + K-) 14.6 1^ 

B* -> J/ipK*, J/ip -> cj)K + K- 9-0 ±i"S 

B± J/ipK±, J/4> -f 2(iT+f£:-) 11.0 ±^*| 
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2.6 io.9 ± 0-3 
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FIG. 3: (a) 2{K + K~) and (b) (j>K + K~ invariant mass spec- 
tra in the r\ c and J/ip regions. The open histograms cor- 
respond to events from the B signal region, and the shaded 
histograms correspond to events from the Mbc-A_E sidebands. 



TABLE II: Measured branching fractions of secondary char- 
monium decays and the world averages Q. The branching 
fractions for modes with K + K~ pairs include contributions 
from -> K+K-. 

Decay mode B (this work) B (PPG) 



rj c — > (p(p 
Vc -> <f>K+K- 
r, c 2(K+K~) 
J/ip -f 4>K + K~ 
J/ip -f 2(K+K~ 



(1.8^I±0.7) 
±1.1) 



+0.9 

0.8 



x 10 
x 10 



(2.9 

(1.4 ±g;| ± 0.6) x 10" a 

(2.4 to.g ± 0.3) x 10~ 3 

(1.4 to B 4 ±0.2) x 10~ 3 



(7.1 ± 2.8) x 10~ 3 

(2.1 ± 1.2) % 
(7.4 ± 1.1) x 10~ 4 
(7.0 ±3.0) x 10~ 4 



Using the known branching fractions B(B ± — > 
J/tpK*) = (1.01 ±0.05) x 10" 3 || and B{B ± -► rjcK^, 
we obtain the secondary branching fractions for J/ip and 
77c decays to 2(K+K-) and <pK + R- listed in TableUD 

Our measured branching fractions for r\ c — > 00 and 
ry c — > 2(i^ + AT _ ) are smaller than those of previous ex- 
periments ||, while those for J/ip decays are consistent. 
The decay r\ c — > 2{K + K~) proceeds dominantly through 
i] c — > <pK + K~ with — > K + K~ . This is the first mea- 
surement of ry c — ► cj)K + K~ . The decay of 77 c — * 00 with 
— * K + K~ makes up approximately 1/3 of the branch- 
ing fraction of r\ c — > <pK + K~ . 

In summary, we have observed the charmless three- 
body decay B — > <p>(pK, which is the first example 
of a b — > sssss transition. The branching fraction 



B(B ± -> 00if ± ) = (2.6+o;g ± 0.3) x 10~ 6 for < 
2.85 GeV/c 2 , is measured with significances of 5.1 a. No 
signal is observed for the decay B — > fj (2220)K with 
/,/(2220) — > 00. The corresponding upper limit at 90% 
C.L. is 5(5* -> /./(2220)/i: ± ) x S(/j(2220) -> 00) < 
1.2 x 10~ 6 . We have also observed significant signals 
for B ± — ► r/cK^ with ?y c — > 00, with ?7 C — > (pK + K~, 
and with ?y c — > 2(K + K~), as well as a signal for B ± — > 
J/ipK^ with J/0 — > cpK + K~ . We report the first mea- 
surement of ?7 C — > (pK + K~ with a branching fraction of 
B(j] c -> (f)K + K~) = (2.9±^;1 ± 1.1) x lO" 3 . Our mea- 
sured branching fractions for ?7 C — > 00 and 2(_?T + i4'~) are 
smaller than those of previous experiments. 
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